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QQ . We study the spectrum of scalar mesons in the holographic dual of walking 



technicolor, obtained by embedding D7-D7 probe branes in a certain type IIB back- 
ground. The scalar mesons arise from fluctuations of the probe techniflavour branes 
and complement the (axial-)vector meson spectra that we investigated in earlier 
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CN ' work. By explicitly finding the spectrum of scalar masses, we show that the non- 

supersymmetric D7-D7 embedding is stable with respect to such fluctuations. In- 
terestingly, it turns out that the mass splitting between the scalar and vector meson 
spectra is of subleading order in a small parameter expansion. It is noteworthy 
that this near-degeneracy may not be entirely due to a small amount of supersym- 
metry breaking and thus could indicate the presence of some other (approximate) 
symmetry in the problem. 
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1 Introduction 

Gauge symmetries play a crucial role in modern particle physics. However, they also imply 
massless gauge bosons. To reconcile gauge invariance with massive gauge interaction 
mediators, one invokes spontaneous symmetry breaking. This is precisely the case with 
the electroweak interactions of the Standard Model of particle physics. Unfortunately 
though, the origin of the symmetry breaking has not yet been established experimentally. 
A leading candidate for the underlying mechanism is the condensation of a weakly coupled 
fundamental scalar field, whose manifestation in nature would be the existence of the 
Higgs boson. There has been a lot of excitement recently about a possible discovery of a 
spin-zero particle at the LHC, which might be the Higgs. However, an actual confirmation 
requires more data. Furthermore, there are other possible explanations for the existence of 
a (fundamental or composite) scalar with mass in the relevant energy range. Regardless of 
the experimental status of the Higgs boson, however, the existence of a fundamental scalar 
would lead to well-know theoretical problems. This motivates the search for alternative 
explanations of the origin of mass in the Standard Model. 



An appealing alternative relies on dynamical effects as the cause of symmetry break- 
ing. In this approach, the role of the Higgs is played by a composite scalar, formed as a 
result of strongly coupled gauge dynamics. Theories of that kind are known as technicolor 
theories [I]. In fact, to be phenomenologically viable, they have to be characterized by a 
gauge coupling that varies very slowly in an intermediate energy range. This phenomenon 
is called walking, in analogy with the more familiar running of couplings. Walking tech- 
nicolor was first proposed in [2]. There has been a lot of work on this subject over the 
years. However, a major hindrance for progress has been the need to study strongly cou- 
pled gauge theories. This is clearly, beyond standard perturbative QFT methods. As a 
result, the computation of electroweak observables from models of dynamical electroweak 
symmetry breaking has long been a challenge. 

In recent years, a new theoretical tool was developed, that addresses precisely the 
investigation of the non-perturbative regime of gauge theories. This is the gauge/gravity 
duality, which allows one to study the strongly coupled regime of a gauge theory via a 
weakly coupled (called dual) gravitational background in a different number of dimensions. 
Thus non-perturbative QFT problems are mapped to (almost) classical gravitational com- 
putations. This powerful method has already produced interesting insights into a variety 
of phenomena, ranging from hydrodynamics [3] to superconductivity [1] . 

Of, perhaps, greatest interest for particle physics is the gravity dual, studied by Sakai 
and Sugimoto |3]. This is a dual of a gauge theory that captures the characteristic 
features of QCD, including chiral symmetry breaking. It arises from a certain D8-D8 
brane configuration in a background sourced by a stack of D4 branes in type IIA string 
theory. The Sakai-Sugimoto model, as well as a wide variety of other similar D-brane 
configurations in both type IIA and type IIB, have been easily adapted to study QCD- 
like technicolor [B]jj In these models one can compute explicitly the S parameter, which 
is an important electroweak observable [9]. However, as is to be expected, the answer is 
not compatible with the known experimental bounds. 

In [TO|[TT|. we studied the gravity dual of a model of walking technicolor. This model 
has two key ingredients. The first is a certain type IIB background [I2] sourced by D5 
branes wrapped on an S"^, whose dual is a walking gauge theory. The second ingredient, 
which introduces the techniflavor degrees of freedom, is a particular D7-D7 probe-branes 
embedding [10] in the background of [12]. In [11], we showed that the S parameter of this 



"'Wc should also note the existence of a huge literature, loosely inspired by AdS/CFT, on bottom-up 
holographic technicolor models [7] , which do not arise from any D-brane construction. 



model can be small enough to be compatible with experiment |j This demonstrates, by 
means of an example, that walking technicolor models can, indeed, satisfy the S parameter 
test for phenomenological viability. 

Here we continue the study of the model of [TOl [TT] . In addition to the spectra of 
vector and axial- vector mesons, that were investigated in [11], there should also be scalar 
mesons. The latter arise from the fluctuations of the technifiavor D7-D7 probe branes. 
As there are two directions transverse to the probe worldvolume, there are two kinds of 
scalar mesons. Since the D7-D7 embedding is not supersymmetric, it is an open question 
whether the scalar spectrum will contain any tachyonic modes or not. In this paper we 
show that the embedding is in fact stable with respect to such fluctuations. Namely, we 
find the mass spectrum of scalar mesons explicitly and thus demonstrate that all mass- 
squareds are positive. Interestingly, it turns out that the scalar meson spectrum is only 
slightly shifted compared to the spectrum of vector mesons|j This is not a priori clear, as 
examplified by the models of [131 E] in which the scalar and vector meson spectra differ 
by order one terms. The near-degeneracy in our case could indicate that the amount 
of supersymmetry breaking, induced by our D7-D7 embedding, is small. Recall that a 
massive A/" = 1 vector multiplet has the same field content as a massless vector multiplet 
and a chiral multiplet. Thus, the massive vector and scalar states belong to the same 
long A/" = 1 representation. However, in the present case, the two kinds of scalar mesons 
turn out to have slightly different mass shifts compared to the vector meson spectrum. 
This suggests that there may be some additional mechanism at work, which would be 
very interesting to understand in the future. 

We should point out that the scalar spectrum investigated here is not directly relevant 
for settling the interesting open question of whether there is a technidilaton. The latter is 
a light composite scalaiu, which is the would be pseudo-Goldstone boson associated with 
the spontaneous breaking of near-conformal invariance of the walking regime. In models, 
in which the breaking of conformal invariance is entirely due to the flavor probes (as in 



^It is, perhaps, worth pointing out that the recent work [8] finds a lower bound on the S parameter in 
a large class of bottom-up holographic technicolor models, which is in conflict with observation. However, 
the considerations in the bottom-up approach rely on certain assumptions, that are not satisfied in the 
actual supergravity duals arising from appropriate D-brane configurations. Thus, the results of [5] do not 

apply to our model. 

^By that we mean that it is shifted only by terms that are sublcading in an expansion in a certain 
small parameter. 

^By light we mean with mass that is (parametrically) much smaller than the masses of the other states 
in the spectrum (in particular, the technirho meson). 



[13], for example), the technidilaton (if present) should appear in the spectrum of scalar 
mesons arising from the fluctuations of those probes. Such mesons are exactly the kind 
of scalar states studied here. However, in our model the technicolor background already 
encodes conformal symmetry breaking. Hence, in our case the technidilaton should arise 
predominantly from the technicolor, not the technifiavor, sectorjfl Thus the question of 
whether our model has a technidilaton is beyond the scope of the present study. We will 
address this problem in another publication. 

The present paper is organized as follows. In Section 2, we briefly recall the basic 
ingredients of the walking technicolor model of [lOl [11]. We also derive the action for the 
two kinds of fluctuations of the flavour probe branes, that give rise to the scalar meson 
spectrum. In Sections 3 and 4 we study the mass spectra of each of the two kinds of 
fluctuations in detail. More precisely, we show that to leading order in a small parameter 
expansion they are the same as the vector meson spectrum of [11]. Furthermore, we 
compute the first subleading correction, that gives the mass shift between the scalar and 
vector boson spectra. Finally, Appendices A and B contain some technical details of the 
computations in the main text. 

2 Action for scalar mesons 

In this section we review the basic ingredients of the technicolor model that we will be 
studying. More precisely, we recall the background that encodes the technicolor degrees of 
freedom, as well as the D7-D7 probe-branes embedding that introduces the techniflavour 
ones. Furthermore, we derive the action for the fluctuations of the D7-D7 embedding, 
which give rise to the spectrum of scalar mesons in this model. 

2.1 Background and flavour probes 

The background we will consider is due to a stack of D5 branes wrapped on an S'^. It 
consists of nontrivial ten-dimensional metric, RR 3-form flux and string dilaton. Further- 
more, it is a solution of the type IIB equations of motion with A/" = 1 supersymmetry. In 
[T2] . this solution was found as an expansion in a small parameter. For our purposes, it 
is enough to consider only the leading order contributions to the background fields. Note 



*It is rather likely that this state will turn out to be related to the light scalar mode found in jl5) . 



that, to this order, the string dilaton is constant. The leading order metric is |12j : 
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The constants c and a above are parameters of the solution and, as in [T2], we have set 
a' = 1 = gs- For an explicit expression for the RR flux, see [12] jj 

This background is dual to a color gauge theory with a coupling constant, which 
exhibits walking behavior. Namely, in a certain intermediate energy range, the gauge 
coupling is almost constant as a function of the energy. In this walking region, p is always 
of order 1 or larger [12]. Hence, coth(2p) ^ 1 while ^^^^(.^^-^ is negligible. In addition, the 
walking region is characterized by [T2] : 



cosh(2p) 



(3 = sin^ a « 1 . 

As a result, to leading order in small /3, one can use the approximations: 
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Therefore, the metric (12. ip simplifies to: 
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To add flavor degrees of freedom, one can introduce D7 probe branes in the above 
background. In fact, as shown in [10], there is an embedding of D7 and anti-D7 branes 



^Note also that the background described above is the leading order solution in an expansion in 1/c 
with c ;» 1. 



such that the D7 and D7 merge at some finite value, po, of the radial coordinate. This 
realizes geometrically chiral symmetry breaking and thus enables us to construct a model 
of walking technicolor. In [TT], we studied the spectra of vector and axial- vector mesons 
in this model. Here we will investigate the spectrum of scalar mesons. The latter arise 
from the fiuctuations of the D7-D7 embedding. So let us, first, describe this embedding 
in more detail. 

The eight worldvolume directions of the probe branes are x'^, p, ip, 9 and (p. The two 
transverse space coordinates are functions of p only: 6 = 0{p), ip = ip{p). The equations 
of motion for 6{p) and (p{p), that arise from the DBI action for the D7 branes, can be 
solved by the following classical configuration fiU\ ITT] : 

±2\^e^P° ^ ^ , (2.8) 

V 3 (1 _ e4(Po-p))i/2 ' ^ ^ 
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dp 
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For future use, le 


coordinate: 









z = ±^f^<^^^^^^V^ ^ (2.9) 

which runs over both the D7 and D7 branes and is such that 2 = when p = pq. Note that, 
up to an overall constant, this is the same worldvolume variable as the one in equation 
(4.2) of [TT]. In terms of z, the classical solution for the derivative of (p is: 

.?-|.?--/|^. (2.10) 

Note also that at z = the above expression for pf is regular, unlike p'^J in (12. 8p at 
p = Pq. This is a key reason to introduce the variable z. Namely, since we would like to 
expand the field p in fiuctuations around the configuration pci, an appropriate choice of 
coordinate should be such that pd does not diverge anywhere along its profile. 



2.2 Fluctuations of D7-D7 embedding 

To find the spectrum of scalar mesons in our model, we want to study 4d spacetime- 
dependent fiuctuations of the D7-D7 embedding around the configuration (12. Sp . So we 
promote the embedding functions to fields in spacetime and expand them as: 

e{z,x>') = e,i + 6e{z,x>') , p{z,x'') = pci{z)+5p{z,x'') . (2.11) 



^"^Of course, we can integrate 'fip{p) to obtain i^(p). However, this will not be necessary or useful, since 
the function Lp{p) will appear in the relevant action below only via its derivative. 

6 



To obtain the action for the fluctuations 59 and 5ip, we need to compute 

C%j = -^-detgsd, (2.12) 



where g^d is the induced metric on the D7 worldvolume, and expand it to second order in 
fluctuationscj To do that, we substitute 

de = e^dz + e^dx'' , dip = (p^dz + p^dx^" (2.13) 



in (12.71) with the change of variable ()2.9p performed. Note that 6'^, 9^ and (y?^ are all first 
order in the fluctuations, since Od = const and d^ipd = 0. However, (p^ does contain a 
zeroth order contribution from (pd- The details of the computation of a/— det gsd can be 
found in Appendix A. The resulting Lagrangian is the following: 



C = const^3{l + z2) [4^2 + gin^^ (4 + c^^ <^2] + q'^po^z^ (4 + c^2 + c(^2) . (2.14) 

To study the action for the fluctuations to quadratic order, we need to expand the last 
expression further. First, however, note that all terms under the square root are quadratic 
in the fields and there is no mixing between the 6 and p> contributions to second order; 
in particular sin^6' = sin^(^ + 59) ^ 1 — (59)'^. Hence, the desired action for each of the 
two fields can be obtained by first turning off the fluctuations of the other field and then 
expanding to second order. We will do this in turn in the next two sections. 

Before concluding the present section, though, let us perform a consistency check on 
the Lagrangian (]2.14p . Namely, let us verify that it does have as a solution of its equations 
of motion the same classical configuration as (12. 8p . (I2.10p . Turning off the x^-dependence, 
which only arises at the level of fluctuations, we have from the square root in (I2.14p : 



'3{l + z^)[9^^ + sm^9p>l] +64^0/3^2 ^ (2.15) 

where we have also divided by an overall factor of 2. The variational equations for 9{z) 
and (p{z) are then: 

Equation (I2.17P is solved by ^ = 7r/2. As a result, (I2.16P simplifies to: 



dz \ Lc 

d {{l + z^)9A (l + z2)sin(2^)((^,)2 



^J?,{l + z^)p)l + e^p^|3z^ 



. (2.18) 



'^We have dropped an overall constant, that is due to the constant string dilaton in this background. 



The latter equation is solved by: 

^^(z) = e'P\l'-^==, (2.19) 



3va 



z 



2 



which is the same as fl2.10p . For completeness, let us also note that fl2.19p implies 

ip(z) = e^P'\- arcsinh(^) , (2.20) 

V o 

where we have set the integration constant to zero in order to have 99 = at 2; = (in 
other words, in order to have the D7 and D7 branches merge at 2; = 0). Finally, note 
that the profile fl2.20p agrees with the one determined by (4.12) of [10]. Namely, in our 
approximations the constant B there becomes B = |; in addition, the plus sign there 
should be taken for 2; > while the minus sign for ^ < 0. 



3 Spectrum of 6 bosons 

In this section we begin investigating the spectrum of scalar mesons in our model. The 
latter arise from the fluctuations of the embedding functions 9 and y? around the U-shaped 
D7-D7 embedding, that realizes chiral symmetry breaking. First, we address the spectrum 
of 9 bosons. In the next subsection we derive the relevant equation of motion and show 
that the resulting mass spectrum differs from the one for vector mesons, computed in [TT] , 
only by an order /3 correction. Finally, in Subsection 13.21 we find explicitly the 0{I3) mass 
difference between the spectra of 9 bosons and vector mesons. 

3.1 Equation of motion 

Let us now turn to studying the spectrum of fluctuations of the field 9{z, x^) around 
the configuration (12. 8p . To do that we have to substitute (12. lip in (I2.14p and expand 
to second order in the fluctuations. Actually, as already pointed out, we can obtain the 
relevant 69 Lagrangian by first setting to zero the if fiuctuations and then expanding. We 
also need to substitute (I2.10p . The result is: 



£ = const ^ {12(1 + z^){dJ9f + e^''°(3 [c(l + z^){d^S9f - 4:{S9f] } . (3.1) 
Using that 3^8^59 = mg69, we then find the following field equation: 



i^WW-3(TT.^ 



{69nz) - ^—^{S9y{z) - J' , 59{z) = MH9{z) , (3.2) 
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where ' = dz and M^ = !!VL^ — . ^jgo^ we have suppressed the x^ argument for brevity. 
To bring (13. 2 p to Schrodinger form, we perform the field redefinition: 

^^^'^ = (TT^ • ^'-'^ 

Then the equation for the mass spectrum becomeso 



where 



Mi + z 



2^2 



VAz) 



e{z) = M'e{z), (3.4) 



Before we turn to solving (13. 4p . let us make a few comments on its form. First of all, 
notice that, except for the Vi{z) correction term, this equation is exactly the same as (4.4) 
of [11] O The latter determines the spectrum of vector and axial- vector mesons in this 
model. Therefore, the presence of the term (13. 5p here implies an 0{P) splitting between 
the masses of 6 bosons and vector mesons. Actually, at first sight one might think that 
such a conclusion is premature, as in [H] we were working at small /3 as well and keeping 
only the leading contributions. So it is legitimate to suspect that the potential in (4.4) of 
[11] might have an order /3 correction that we have just not written down. However, we 
will show now that any such correction would be the same for both the vector bosons and 
6 mesons. Hence the correction (13.51) does, in fact, give rise precisely to the mass splitting 
between the 6 and vector meson spectra. 

To understand that, it is most instructive to derive the full field equation for 66, 
without using any of the approximations that led to the simplified walking metric (12. 7p . 
Working with the full metric (12. ip and introducing the same functions a{p) and b{p) as 
in (2.21) of [n], namely 

b{p) = ^[m+gipHY^'g"' (3.6) 



^^Lct us make an important remark. Clearly, the potential term in the Schrodinger equation (|3.4[) 
changes sign as z is varied. However, the presence of a negative potential region does not by itself 
indicate an instability. The configuration we are expanding around would be unstable if the spectrum of 
solutions of p.4p with appropriate boundary conditions were to have a negative mass-squared mode. In 
the rest of this Section, we will show that this is not the case. 

^■^In that regard, note that the parameter A in (4.5) of [11] can be removed by introducing a new 
variable z = z/^/X in (4.4) there. In addition, using that co = c^/P and the expression for A in (4.6) of 
[TTj . one can see that the new variable z is precisely the same as the variable z in (|2.9p here. 
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48 4 48 4 


25Q n (Pjn (PJ coth (.pj , g[p) 2^256 




and r ^'^'^'^>'^yP^ 



128 (/ + g^l) 
we find after a lengtliy calculation: 



Pi'(p)P^(p) coth(2p) (3.7) 

(3.8) 



- d, [bip) d,5e{p) ] - h{p) coth^(2p)b^'(p)]25e(p) = ml a{p) Se{p) . (3.9) 

The above equation of motion differs from the (axial-)vector meson one, eq. (2.26) 
of [TT], only by the presence of the [v^^']^ term. Although it is by no means easier to try 
to solve (13. 9 p instead of (13. 4p . it is instructive to see how (13. 9p reduces to (13. 4p in the 
approximations that give the walking metric (12. 7p . Namely, to leading order in small /3, 
we find that the functions (13. 6 P reduce too 



3l/V/4/33/V° (1 + ^2)3/2 

^^^ '' 24 z 

KpW) = TT^ ^^, (3.10) 



VT 



z 



2 



where we have used (12. 9p . It is then easy to see that (13.90 becomes exactly (13. 2p . after 
the change of variable from p to z. The important point is that the term — 3^+^*2 ■) '^^(-^) 
in (13. 2p . that gives rise to Vi{z), comes entirely from (and in fact, to leading order, is 
the only contribution of) the [v^p]^ term in (13. 9p . which is precisely the extra term not 
present for vector bosons. Hence, whatever other 0{P) contributions might arise from the 
remaining terms in (13. 9p . when subleading orders in a{z) and b{z) are taken into account, 
they will necessarily be the same for both the vector and 6 mesons. Thus the leading 
order mass difference between the two spectra is determined precisely by Vi{z), while any 
subleading correction to it will come from higher order terms in the small f3 expansion of 
the [y^p]^ term. 

3.2 Mass spectrum 

Let us now consider the Schrodinger equation (13. 4p in more detail. Clearly, at leading 
order in small /3, the Vi{z) term in the potential can be neglected and so one has exactly 



^''^To obtain p.lOp . one also needs to use the relation cy^ = -j^, valid to leading order in small /?, that 
was derived in 1111. 
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the same spectrum as for the vector mesons, studied in [TT]. To evaluate the mass spht 
between the two spectra due to Vi{z), we will use perturbation theory. Note also that, 
since we are studying the perturbative spectrum, we will impose the boundary conditions 
that at z = ±z\ the wave function vanishj^i Let us denote by Q^ the zeroth order wave 
function of the n^^ state, i.e. with Vi{z) neglected. Then, according to standard rules, 
the correction to the mass-squared M^ due to the perturbation Vi{z) is given by: 



^\Ml] = —J [Ql{z)fV,{z)dz, (3.1i: 

J-ZA 



where the normalization factor is 



fZA 

N'= I [Ql{z)fdz. (3.12) 



ZK 



From the analysis of Section 4 in [TT], we know the approximate analytical form of 
Q^{z). More precisely, one can divide the interval (0, za) into two regions, in each of 
which the field equation simplifies and can be solved analytically. The matching of the 
two then gives the quantization condition M^ = rn/z^, where r„ is the n**^ root of the 
Bessel function Jq- Up to overall constants, the wave function has the following form in 
each of the two regions: 

no. ^ j <dl:\z) = {I + z'y^\ for smalls, i.e.^e(0,^.) , 

[^6^^(2;) = A/iJo(M„2;) , for large z, i.e. z e {z,,zt^ , 

where 1 « 2, « -zaUj Actually, it is easy to realize that one can approximate the solution 
on the entire < z < za interval (except for a A2; ~ M^ neighborhood oi z = \fi) by the 
following function: 

e°(^) = (l + z2)V^Jo(M„^). (3.14) 



^^Here, as in |11| . we have introduced a physical cut-off za at the upper end of the walking region, so 
that z G (— 2a,2:a). The earlier work |10| discussed holographic renormalization for the background of 
interest. However, in [TT] we pointed out that above the scale za the light spectrum of the model should 
change in order to reflect the degrees of freedom appropriate for extended technicolor. Hence, on physical 
grounds, our current model should not be viewed as valid to arbitrarily high energies. The proper UV 
completion is still an open problem; see j28| . though, for important recent progress in that direction. 

^^For simplicity, we took in p.l3p the form of the solution, which is symmetric under z — >■ —z. The 
antisymmetric solution, that now gives pseudoscalar mesons, is a bit more involved technically as can 
be seen from the axial- vector considerations in [11[ . However, the final outcome is the same, except for 
substituting r„ with the quantity [in defined in (4.16) of [TT] . 
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Clearly, reduces to the first line of (13.131) at small z, while it acquires the form on the 
second line of fl3.13p at large z. To see by how much this approximation fails to be an 
exact solution, let us substitute = in (13.41) . with Vi omitted. We obtain: 



Ql"{z) 



^ -2 ,,2 



A{l + z 



2^2 



Qliz) 



MnJljMnZ) 
■ z(l + ^2)3/4 



(3.15) 



where the nonvanishing of the right hand side is exactly the measure of the deviation of 
from being an exact solution. It is easy to see that, aside from a small neighborhood 
of a finite number of points, this right hand side is of 0{z~^) compared to the individual 
terms on the left hand side on the whole interval < 2; < za- 

In evaluating the integral in (I3.12p . it is easy to realize that the leading contribution 
in the small /3 expansion will come entirely from the large z region. So, to leading order, 
we have: 



N' 



2 A r^K 



[e:{z)Ydz 



ZK 



z[Jo{Mr, z)Y dz = zi [Ji(M„2;a)]' • 



(3.16) 



On the other hand, the integral in (13. lip can be estimated by using the approximate 
solution 0°(2;). To leading order in z\, one obtains: 



ZA n 

[Ql{z)rV,{z)dz = —Pe^'"' 



arcsinh(2A) - -r^ sF,{l, 1, -; 2, 2, 2, 2; -r^) + dz], 



(3.17) 

where 3-^4(1, 1, |; 2, 2, 2, 2; — r^) is the hypergeometric function. The details of this com- 
putation are somewhat subtle and we have relegated them to Appendix B. Combining 
the above results, we find to leading order in z\: 



S[M'n 



13 e 



4po 



3 r„Vi2(r„; 



Mf 



arcsinh(2;A) - -fi 3^4(1, 1,^2, 2, 2, 2; -r^) + 0{zl' 



(3.18) 

where we have used za = ^n/Mn. Now, recalling that M^ = !!!:ii^ — and also using that 
for large zx the function arcsinh(2A) = ln(2;A + a/1 + z\) simplifies to ln(22A)5 "we can 
write the correction to the boson mass spectrum, i.e. 5{mP^Y = {m^Y — (m^)^, as: 



e\2 



KO 



m^ 



2f3e 



4po 



3[r„Ji(r„)]2 



log(2zA) - -In 3^4 



r! 3^4(1, 1,^:2, 2, 2, 2: -r!^ 



^) -^) Q ) ^^ ") ") ") 



0{z-, 



where the unperturbed (i.e., with Vi{z) neglected) mass of the n*^ state is [TT 

12 



m„ 



Co 



:^-i/2g-4p^ 2^2 



n^^-'-KK 



(3.19) 



(3.20) 
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with Co = ^ to leading order. 

As can be seen from f l3.19p . the mass spht between the n*'^ vector and 6 boson states 
is suppressed by a power of /3 compared to the leading contribution. It would be very 
interesting to understand the origin of this split, or rather the reason for its suppression. 
We leave this for future work. 

4 Spectrum of Lp bosons 

In this section we will study the fluctuations of the field ip{z,x^) around the classical 
solution (12. Sp . ( I2.10p . As before, to obtain the Lagrangian of interest we have to substitute 
(12. lip in (I2.14P and expand to second order in 5ip{z, x^) and 59{z, x^). However, as already 
pointed out, to quadratic order there is no mixing between the 9 and c/? fluctuations. So 
we can safely set 9 = n/2 from the start. Hence, the Lagrangian (I2.14p acquires the form: 

C = constsj 12{1 + ^2) ((^'^^ + 5^if + e^po/? ^2 [4 + ^ {d^5^f] , (4.1) 

where again ' = d^- Expanding to second order and substituting (I2.10p . we have: 



Se-^"" ^' ..„n2 , o,A7x„. , ^"""^^ 



2 



where for convenience we have dropped the overall constant in front of the square root 



in (14. ip : also, we have omitted the additive term 2e'^^°^/^^/l + z^ , since clearly it drops 
out of the field equation. Now, varying (14. 2 p and using d^d^5ip = m'^dip, we obtain the 
following equation of motion: 

- V(^) - 4^ ^^'(") = K ^^(^) ' (4-3) 

where M^ = im? c^e^''^'' / 12 and we have suppressed the argument x^ for brevity. 

In line with the previous section (as well as our considerations in [H]), it would seem 

that the strategy to study (14. 3 p invloves first transforming to Schrodinger form. For 

technical reasons, however, we will use another approach. Despite that, let us briefly 

discuss for completeness the Schrodinger form of equation (14. Sp . It is achieved via the 

field redefinition: 

5^{z) = ^ ^'> ^z) , (4.4) 

z 



which leads to: 
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This is rather similar, although not quite the same as the equation for the vector mesons. 
Recall that the latter is 

-^'\z)-^^^^^{z) = Mlnz), (4.6) 

as explained in the previous section. There we found that the difference between the 9 
and vector meson cases was given by the order (5 correction to the potential Vi{z) in (13. 5p . 
On the other hand, (14. 5 p and (14. 6 p differ by a potential difference ^V = ■7^^2)2- At large 
z, one has AV" ~ ^, which is negligible compared to the large z behaviour of each of the 
two potentials in (14. 5 p and (14. 6p . that is ^. However, at small z the difference is of the 
same order as the potentials themselves. So we cannot view AV as a small perturbation 
compared to the vector meson case. Nevertheless, one may expect that the mass difference 
between the two cases is small due to the same large z asymptotics they exhibit. 

Since we cannot solve exactly either of (14.50 and (14. 6p . it turns out that the most 
efficient way to compare their mass spectra is to take a step back and to not go to 
Schrodinger form. That is, for the ip bosons we will consider (14.30 : for future use let us 
also introduce the notation: 

H - ^' ^ + ^' ^ (A7) 

^ dz^ z{l + z^)dz ^ ' 

for the differential operator on its left-hand side. Clearly, H^p can be viewed as a Hamil- 
tonian operator with eigenfunctions 5ip{z) and eigenvalues M^ . The analogue of (14. 3 p for 
the vector bosons is 

-<(z)--^^'(z) = M^^(z), (4.8) 

where 

V'(-) = (T^^ (4-9) 

with the same '^{z) as in (14. 6p . Note that in [H] we did not write down equation (14. 8 p 
explicitly. Instead, we obtained directly (14. 6p . because we combined the change of variables 
from p to 2 with the redefinition (14. 9p . Indeed, the factor ., .^. .,^/^ in (4.3) of [11] accounts 
precisely for the factor L^^^ in (14. 9 p here. As a last preparation, let us also introduce 
the notation: 

*- = -dI5 - TTT-Tz '^■^°' 

for the differential operator on the left-hand side of (14. 8p . 

Now, our goal will be to compare the Hamiltonians H^p and Hy- Since, as already 
pointed out above, the mass difference is expected to be small, we can view Aif = H^—Hy 
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as a small perturbation. Of course, once we obtain the answer, it will be easy to verify 
that this assumption is satisfied. Before proceeding further, it is very useful to notice the 
following relation: 

H^-Hv = 2{H^ - Ho) , (4.11) 

where the operator Hq is 

dz"^ z dz 
The differential equation Hqx{z) = M^x(z), namely 

-X"i^)--X{^) = M\{Z), (4.13) 

z 

is in fact precisely of the form that both (14. 3 p and (14.81) acquire at large z. Unlike those 
equations, though, (14.13^ can be solved easily. The regular solutions, that vanish at 
z = ±2a, are given by: 

Xn(^) = Jo(^), (4.14) 

where r„ is again the n^^ root of the Bessel function Jq. The key use of (14. lip is therefore 
that, since we know the eigenfunctions of Hq, we can compute the mass shift, due to AH, 
by using standard first order perturbation theory: 

AM2 = -^- {Xn\AHo\Xn) , (4.15) 

\Xn\Xn) 

where AHq = H^p — Hq. And, of course, since (I4.13P is the same as the large z asymptotic 
form of (14.30 . the perturbation AHq is expected to be small just as AH = H^ — Hy 
is. Before turning to the computation of AM^ according to (I4.15p . we also need to note 
the following. Hermiticity of the Hamiltonian Hq requires the use of the measure \z\dz. 
Indeed, one can easily show that 

Xm{z)Xn{z)\z\dz r^ 5mn (4.16) 

-ZA 

and 

ZA j-ZK 

Xm{z)HoXn{z)\z\dz = / Xn{z) Hq Xmiz) \z\ dz ~ 5mn ■ (4.17) 

Zjy J-ZA 

Now we are finally ready to compute AM^ = [M^]^ — [M^]^. Substituting the explicit 
form of Xn{z) in (I4.15p . we find: 

4 r^^ 1 d 

[M^f-[M^f = -—j^ zJ,{M^z)-^-—^^-J,{M^z)dz, (4.18) 
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where 



m 



2A 



z[jQ{M^z)fdz = zl[Ji{rn)f 



(4.19) 



On the last line we have used Jo(M„za) = 0. After substituting x = M^z , we obtain: 



yi2 



[Mn-MY = -—, 



m 



Ml 



d 



M2 + x2 dx 



[Jo(x)]^(ix 



m 



Ml[\og{l/Mn) + Cn + 0{Ml log(MO)] 



where the constants c„ are determined by: 



1 /■'■" 1 



K • ■■' 2 



X 



r^{x) + --i 



dx . 



(4.20) 



(4.21) 



Equivalently, we can rewrite O4.20p as: 



A(m^) 



f\2 



(m, 



ip\2 _ 



m. 



V\2 



m^ 



:[\og{zJrn) + c„ + 0{zl^ \og{zf,) )] , (4.22) 



where again rn^ is given by (I3.20p . For example, for the lowest lying state ci ~ 0.619. 
Thus, for n = 1 we have: 



?Tif )^ — [m 



V\2 



3.86 



m^ 



[log(^A) + 0.887 + 0(^XMog(^A) )] 



(4.23) 



It is interesting to compare the mass differences between the two kinds of scalar bosons, 
9 and v?, and the vector mesons. First, from (I3.19P and 04.22p one can see that A(?7i^)^ 
and A{m'-^Y have opposite signs. More precisely, the 9 bosons are slightly lighter, while 
the ip mesons are slightly heavier than the vector bosons. Despite that, the magnitudes 
of A(?7i^)^ and A(?72^)^ are almost comparable. This is because /3 



^A ' 



from za = z{pf 



which follows 
z^p'^ as well as the relation e^^^ ~ /3~^/^ that was elaborated on in 
[TT]. However, to be more precise, we should write e'^''^ = o"/3~^, where a is a number of 
order 10"^ or 10^^ \]A\. Therefore, the ratio A(?72^)^/A(?7i5^)^ ~ a < 1 and thus the if 
boson shift is somewhat greater than the 9 boson one. 



5 Discussion 



We investigated the spectrum of scalar mesons in the gravity dual of a walking techni- 
color model, obtained by embedding D7-D7 techniflavor probes [ini [H] in the technicolor 
background of [12]. We showed that all mass-squareds are positive and thus there is no 
tachyonic mode with respect to fluctuations of the probe branes embedding. Furthermore, 
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we found that both types of scalar mesons, arising from those fluctuations, have spectra 
that are nearly degenerate with the vector meson spectrum, studied in [11] O Namely, the 
mass differences are of higher order, in an expansion in the small parameter /3, compared 
to the leading contributions to the masses. It would be very interesting to understand the 
underlying reason for this near-degeneracy. Likely, it will mostly turn out to be due to a 
small amount of supersymmetry breaking, as the massive vector and scalar states should 
belong to the same long A^ = 1 supermultiplet^j However, the difference between the 
mass shifts of the two kinds of scalars may indicate a role for an additional factor. We 
leave investigation of this question for future work. 

One should keep in mind that the scalar spectrum studied here does not address the 
question of whether our walking technicolor model has a dilaton or not. The latter is 
also a scalar field. However, it is a mode that, for us, should arise predominantly from 
fluctuations of the technicolor background (see the discussion in the introduction). In 
particular, its distinguishing feature in the effective Lagrangian is a certain coupling to 
the square of the technicolor field strengthO The phenomenological importance of this 
field is great, since it is a (potentially) light scalar boson rather similar to the Higgs. To 
differentiate between the two, one would have to compute their couplings to other fields; 
for such a study see [H] . There is a lot of literature for [191 [201 |2ll 122] or against [231 [Ml [H] 
the existence of a dilaton in walking technicoloro The arguments for are inspired by the 
approximate conformal symmetry of the walking region, whereas the arguments against 
rely on the fact that the theory is never exactly conformal. We hope to investigate the 
issue of the presence of a technicolor dilaton in our model in a future publication. 

Other important open problems are the computation of the T and U electroweak 
observables, as well as the anomalous dimension of the technifermion condensate. Finally, 
there are various interesting possibilities for extending our model. An important issue, 
equally relevant for all other D-brane constructions of technicolor models [6], is to go 
beyond the flavor probe approximation. This would enable one to study gravity duals of 



^^Interestingly, this is in line (albeit, more pronounced in our case) with the expectations, based on 
the use of the Schwinger-Dyson and Bethe-Salpeter equations, that in IR-fixed-point walking theories the 
masses of scalar and vector mesons at a given level are numerically close to each other [TB] . 

^^Such an expectation is supported by the results of ^7\, where the flavor probe embedding is super- 
symmetric and the spectra of vector and scalar mesons are exactly degenerate. We thank C. Nunez for 
pointing this out. 

^^Note that [15] found evidence for a light scalar in the spectrum of fluctuations of the color background 
of [12]. However, at present, it is unclear whether this state can be identified with the technidilaton. 

^°It is also worth mentioning the pioneering work of [25] on Goldstone bosons for broken conformal 
invariance, albeit in the case of weakly-coupled gauge theories. 
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models with comparable numbers of colors and flavors, like the walking technicolor models 
that rely on the existence of a Banks-Zaks flxed point for the gauge coupling in the infrared 
region [26]. Another interesting direction, relevant for understanding the generation of 
quark and lepton masses, is to try to build a gravity dual of extended technicolor [27] . In 
that regard, the recent work [28] might provide valuable insights. 
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A Determinant of induced metric 

Here we compute the determinant of the fluctuated metric, induced on the D7-D7 world- 
volume. 

Performing the change of variable (12. 9p , the metric (12. 7p acquires the form: 



/7s2 _ A 



cP[ 



r]^,dx^dx'' + ^ (4 {d,pf dz" + (^3 + Co-if) + ^ [d^l + d^^) 



(A.i: 



where recall that ws = cos 9 dip and dfig = d9 + sin 9dip , also, now we have 



P[ = l(3e''P^^(z'' + l). (A.2) 
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Substituting (12.131) in (]A.1|) . we find the following non- vanishing metric components: 



cP' 

'n^lu + -^{cos^ 0)ipf,ip^ + '-O^Ou + ^(sin^6')v9^v9j. 



9iiz = A 
9zz = A 

9zLb-i 






cP' 



'-^md,pf + cos'e^l) + '-{9l + sm'e^ 



A——- coso' (fz 



9. 



UJ3UJS 



A- 



cPi 



9lJ.U]3 

c 



A—— cosO ip^ 



fl'tJilii ^ 



9i: 



A- 



(A.3) 



4 ' ^'^^'^^ 4 

Hence, we have that det gsd = ^ (f) dett^ed, where g^d is the 6d metric along the yU, p 
and 0)3 directions. 

Now, we can use that 

det g^d = 9cjzGji det G^d , (A.4) 

where G^d is a 5d metric with the following componentscj 

9Mu)39NLb-i 



Gmn — 9mn — 



9i 



with M, A^ running over /i, z 



U3UJ3 



Hence 



A. 



detgsd = A''(^j) I^detGsrf 
where, more explicitly, the components of G^d are: 

cP' 

G^u = 9iiu - A—^ cos^6' (/5^v9, 

o 

cP' 

G^,z = 9i^z - A—^ cos'^9(ff,Lpz 

o 

cP' 

Gzz = 9zz - A-^ cos'^Oipl . 

Substituting g^^, g^z and gzz from ( 1A.3I) . we find: 



(A.5) 



(A.6) 



(A.7) 



G 



flU 



A 



VtMu + -. {d^fiu + sin^6'(^^(/)^) 



G^z = A-{e^ez + sm^eip^ipz) 



G, 



A 



cP[ 



^2 , C 



2 {dzpf + \{el + su.^9^l) 



(A.8) 



^^This is, clearly, a special case of the general relation det = det A x det (Z? — CA ^i?) 
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Now we can single out the z coordinate and write det G^d = Gzz det G^d with G^d 
determined similarly to G^d in (lA.Sp . Therefore, we have: 



/ c \ 2 cP' 
det^8<i = ^^(^j -^CdetG4d: 

where Gm is a 4d metric with the following components: 



(A.9) 



'^flU ^ ^v 



Gzz 



(A.IO) 



In order to compute the contribution of the second term in ( lA.lOp . note that we want to 
expand to second order in fluctuations. Given that the classical solution corresponds to 
= const, all derivatives 6'^, 6z are first order in the fluctuations. However, the classical 
solution for if depends on z. So, while v?^ is first order, the derivative ipz has a zeroth 
order (i.e., classical) part. So, expanding (lA.lOp to quadratic order in the fluctuations, 
we have: 



'^ flU ^IJ,U ^ 



vl 



? {dzpf + f V^^ 



Vt^^y 



(A.ll) 



2 y-^zr, I 4Y-Z 

where in the coefficient of (p^i^y we have taken the classical part of (pz and substituted 
the classical solution siv?6 = 1. Hence, to quadratic order, the determinant of G4d is: 



- det Gm = A^ 



1 ^ ^ r.A-/Q ft ^ g-Pi i^zp) 

'^4^ '^''^2[2Padzpf^,l\ 



^^"Vi.^u 



(A.12) 



Combining this with (IA.9P and (lA.Sp . we finally find: 



det^8d = -^*^( 



c\^ cP[ 



'-^{dzpf + 1(91 + sin^e^l) 



{^ + l^l + Qvl), (A.13) 



where 



Q 



cP[ {dzpf 



(A.14) 



2[2P[{dzpf + ^l\' 
Therefore, substituting (lA.2p and the explicit form of dzp from (12. 9p . we obtain the fol- 



lowing expression for the Lagrangian C = — \/ — det g^d '■ 



C 



96 



3(1 + z2) [4 ^2 + sin2^ ^4 ^ ^ ^y ^2] + ^Ap^p ^^A + cOl + c ^l) . 

(A.15) 
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B Estimating boson mass shift 

In this Appendix we derive the result f l3.17p for the integral in the numerator of (13. lip . 
We use the approximate wave function f l3.14p . whose failure to be an exact solution is of 
order z^^. Hence we have: 



where 



[Ql{z)]'V,{z)dz = --/5e^^«/, (B.i: 



rzA 1 

/ -j===[MM^z)]Uz . (B.2) 

Jo Vl + -2 



To evaluate /, let us rewrite it as: 



r^r^+ r ^^ {[MMnz)]' - 1} dz . (B.3) 

Jo V 1 + 2"^ Jo Vl + z'^ 



The first integral is easy to compute: 



rzA I 
h = —^^=^= dz = arcsinh(zA) . (B.4) 

Jo V 1 + z^ 



To estimate the second integral in (IB.3p . let us first substitute x = MnZ: 

dx 



'0 VM2 + x2 

Although this cannot be solved analytically, we will extract now its leading contribution 
in a small parameter expansion. For that purpose, recall that M^ = j^e^''° f3^^'^m'^ with 
the unperturbed (i.e., without Vi) mass-squared being m,^ ~ (9(/3^/^) [11]. Hence M^ is 
of order /3 and so can be viewed as a small parameter. Then, at zeroth order in /3, we 
have: 

/ f {['^o(x)F-l} = -4r^3F4(l,l,|;2,2,2,2;-0, (B.6) 

which notably is independent of z\. To estimate the error in omitting M„, let us rewrite 
I2 identically as: 



/o VM2 + X2 i^ 2 J 2^0 ^/W+^^ 

This is useful for the following reason. Let us consider the first integral in (IB.7P and recall 
that for small argument [Jo(x)]^ f« 1 — |x^ + 0{x^). Then that integral does not have a 
logarithmic divergence in the limit M^ — )■ and, furthermore, its linear term in a small 
M^ expansion is finite. Hence, we can obtain its contribution to 0{M^) by expanding the 
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integrand in small M^. On the other hand, the second integral in (IB. 71) can be performed 
exactly. Combining these two contributions and substituting M„ = r„/zA, we find: 

h = -\rl 3^4(1, 1, 1 2, 2, 2, 2; -r^) + ^ log(^A) + 0(^^2) . (B.8) 

The sum I = Ii + I2 then gives precisely (13.171) . 
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